/* This file is part of the Palabos library.
 * Copyright (C) 2009 Jonas Latt
 * E-mail contact: jonas@lbmethod.org
 * The most recent release of Palabos can be downloaded at 
 * <http://www.lbmethod.org/palabos/>
 *
 * The library Palabos is free software: you can redistribute it and/or
 * modify it under the terms of the GNU General Public License as
 * published by the Free Software Foundation, either version 3 of the
 * License, or (at your option) any later version.
 *
 * The library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

/* Main author: Orestis Malaspinas
 */

/** \file
 * Descriptor for all types of 2D and 3D lattices. In principle, thanks
 * to the fact that the OpenLB code is generic, it is sufficient to
 * write a new descriptor when a new type of lattice is to be used.
 *  -- generic code
 */
#ifndef ADVECTION_DIFFUSION_LATTICES_HH
#define ADVECTION_DIFFUSION_LATTICES_HH

#include "latticeBoltzmann/advectionDiffusionLattices.h"

namespace plb {

namespace descriptors {

    // AdvectionDiffusion D2Q5 //////////////////////////////////////////////

    template<typename T>
    const T D2Q5Constants<T>::invD = (T)1 / (T) d;

    template<typename T>
    const int D2Q5Constants<T>::vicinity = 1;

    template<typename T>
    const int D2Q5Constants<T>::c
            [D2Q5Constants<T>::q][D2Q5Constants<T>::d] =
        {
            { 0, 0},
            {-1, 0}, {0, -1}, {1,0}, { 0,1}
        };

    template<typename T>
    const int D2Q5Constants<T>::cNormSqr[D2Q5Constants<T>::q] =
        { 0, 1, 1, 1, 1 };

    template<typename T>
    const T D2Q5Constants<T>::t[D2Q5Constants<T>::q] =
        {
            (T)1-(T)2/invCs2, 
            (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2), 
            (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2)
        };

    template<typename T>
    const T D2Q5Constants<T>::cs2 = (T)1 / (T)3;

    template<typename T>
    const T D2Q5Constants<T>::invCs2 = (T)3;

        
    // AdvectionDiffusion D3Q7 ////////////////////////////////////////////////////

    template<typename T>
    const T D3Q7Constants<T>::invD = (T)1 / (T) d;

    template<typename T>
    const int D3Q7Constants<T>::vicinity = 1;


    template<typename T>
    const int D3Q7Constants<T>::c
        [D3Q7Constants<T>::q][D3Q7Constants<T>::d] =
        {
            { 0, 0, 0},

            {-1, 0, 0}, {0,-1, 0},
            { 0, 0,-1}, {1, 0, 0},
            { 0, 1, 0}, {0, 0, 1},
        }; 

    template<typename T>
    const int D3Q7Constants<T>::cNormSqr[D3Q7Constants<T>::q] =
        { 0, 1, 1, 1, 1, 1, 1};

    template<typename T>
    const T D3Q7Constants<T>::t[D3Q7Constants<T>::q] =
        {
            (T)1-(T)3 / invCs2,

            (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2),
            (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2), (T)1/(invCs2*(T)2)
        };
        
    template<typename T>
    const T D3Q7Constants<T>::cs2 = (T)1 / (T)4;

    template<typename T>
    const T D3Q7Constants<T>::invCs2 = (T)4;

}  // namespace descriptors

}  // namespace plb

#endif
